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We find a modified Bloch equation for the electronic magnetic moment when the magnetic 
moment explicitly contains a diamagnetic contribution (a magnetic field induced magnetic moment 
arising from the electronic orbital angular momentum) in addition to the intrinsic magnetic moment 
of the electron. The modified Bloch is coupled to equations of motion for the position and momentum 
operators. In the presence of static and time varying magnetic field components, the magnetic 
moment oscillates out of phase with the magnetic field and power is absorbed by virtue of the 
magnetic field induced magnetic moment, even in the absence of coupling to the environment. We 
explicitly work out the spectrum and absorption for the case of a p state electron. 



I. INTRODUCTION 

The importance of magnetic resonance methods to chemical analysis, solid-state physics, biological structure de- 
termination and medicine is well recognized. In magnetic resonance experiments, a material is placed in dc and 
oscillating magnetic fields and the power absorption is measured. The power absorption exhibits spectra containing 
resonances, from which valuable information about the system is obtained. 

The Bloch equation describing the dynamics of the magnetic moment operator fl is: 

Its solution is discussed in textbook theoretical treatments of magnetic resonance Here, /io — ^ is the Bohr 
magneton, and 7 is a dimensionless constant. SI units are employed throughout. This equation has played a major 
role in understanding magnetic resonance experiments ever since it was first derived by Bloch in 1946 j^]. In the 
usual treatment, coupling to the surrounding, i.e., to additional degrees of freedom of the system (often described as 
a thermodynamic bath), is incorporated into the description by adding phenomenological relaxation (decay) terms 
to the Bloch equation . The additional phenomenological decay terms arise when the bath degrees of freedom 

of the system are adiabatically eliminated from the description. The external magnetic field in magnetic resonance 
experiments has a temporally-constant strong component, and a weaker oscillating component of radio- frequency or 
microwave radiation, which in one common configuration is linearly polarized and perpendicular to the dc component: 

B = Boz + Bi {e"^' + e"*"*) x ; Bi < Bq ■ (2) 

It is a straightforward generalization to treat the case when the magnetic field contains a rotating oscillating field: 
Be = Bqz ^/2Bi[x cosiujt) — ysm[LLit)\. Here we consider linear polarization. For atomic and molecular systems it 
is an excellent approximation to neglect the spatial dependence of the oscillating field because the wavelength is so 
much larger than the size of the atom or molecule. 

The power absorption P of a magnetic moment pL in the external field B is given by: 



P 



dB 



d{ll) 



dt 



(3) 



where (...) stands for quantum expectation value and [...J^ for time averaging. Power can only be absorbed if B and pi 

oscillate out of phase, and fl has a component along The Bloch equation without relaxation terms, Eq. (|^), results 

in no power absorption, since clearly, B ■ {j3 x pj =0. In this treatment, the phase shift between the oscillations of 

B and fl, necessary for power absorption, exists only because of the relaxation terms. Without those, Eq. (^, implies 
that power absorption is impossible. 



1 



This theoretical framework is adequate for the description of both NMR (Nuclear Magnetic Resonance) and ESR 
(Electron Spin resonance). It is generally not correct for EMR (Electron Magnetic Resonance). The Bloch equation 
is not the correct equation of motion for the magnetic moment operator of an electron with non-vanishing orbital 
angular momentum. Orbital angular momentum effects in magnetic resonance phenomena can be important, not 
just in gases, but even in solids, e.g., the transition and lanthanide ions doped in crystals where quenching (the 
modification of the magnetic susceptibility to values close to those produced by electron spin alone due to large 
crystal field splittings) does not occur §. Here we show that, due to diamagnetic effects arising from the electronic 
orbital angular momentum, the Bloch equations are replaced by a set of coupled operator equations for the magnetic 
moment jl (which incorporates a magnetic field induced magnetic moment contribution), the coordinate r and the 
velocity operator Il/m, defined below. The interaction Hamiltonian has non-linear terms arising from a magnetic field 

induced magnetic moment. As a result, the magnetic moment has an oscillatory component along ^ and oscillates 
out of phase with the external field. Hence, power absorption does not vanish, even when coupling to the environment 
is not included. 

In Sec. H we define the magnetic moment operator, including the magnetic field induced magnetic moment contri- 
bution, and in Sec. Ill we derive the Heisenberg equations of motion for a magnetic moment pi in the magnetic field 
B. Three coupled equations are obtained in this way, Eqs. (|Tl|), ([T^), and ( p^ below. Based on these equations, 
Sec. 1^ develops general considerations of power absorption in EMR experiments. Since, unlike the Bloch equation, 
these coupled equations can not be solved in closed form, we solve them in time-dependent perturbation theory, and 
calculate the power absorption using this analysis. In Sec. ^ we focus on the example of an electron bound in a 
p-atomic-orbital, and derive Eq. ( p3| ) for the power absorption to leading order in the fields for this case. Section VI 
presents a summary and conclusion. 



II. THE HAMILTONIAN AND THE MAGNETIC MOMENT OPERATOR 

The Hamiltonian for an electron, with charge q = —e and mass to in a spherically symmetric potential V{r) and 
an external magnetic field B{t), can be written in terms of the kinetic velocity operator H/m ||^: 

H-nr)^^^ (4) 

Il = j>-qA, (5) 
A=]^By.'r. (6) 



For simplicity, we have neglected the electron spin (which can be re-introduced rather easily as shown in Sec. VI) 
and chosen a symmetric gauge. The coordinate operator r, the conjugate momentum operator p, and the velocity 

operator H/to, satisfy the commutation relations: [ri,pj\ — [r.i,Hj] = ifiSij, and [Hi,Hj] = ihqeijkBk, where eijk is 
the anti-symmetric Levi-Chivita tensor. 

The gauge-invariant operator for the orbital magnetic dipole moment of the electron is given by 

^=J_;-'xri. (7) 

The magnetic moment operator defined in Eq. (0) is gauge invariant. This can be seen by considering the gauge 
transformation matrix U = exjp[i^x\ in terms of the gauge function x{r,t). Under the gauge transformation, the 

vector potential and the magnetic moment are transformed as follows: A ^ A' = A + Vx, and fl ^ jl ^ fx 
(p'-gi')]t/t=^. 

The Hamiltonian of Eq. (Q) can be written in terms of the magnetic moment in the fohowing form: 

"2 2 

H = V{r) + ^-p.B + -^\{B-f)^ - B^rA . (8) 
2to 8to L J 

Neglecting all terms of order B'^ gives: H' — V{r) + 4^ — fJ-' ■ B and fi' = f x p, as used in standard textbook 
treatments of magnetic resonance. Here, we consider the effect of not neglecting these second order terms, i.e., using 

H (not H') with the full /2 (not /T'). 
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III. HEISENBERG EQUATIONS 



The magnetic moment fl satisfies the Heisenberg equation of motion: 



dfi 
It 



djl 
dt 



gives, 



dpi 
'dt 



2m at 



4m 



4m 




Straight-forward algebra 



,d B 
' dt 



(9) 



\(r-^p)B-(r-B)p + 
2mV {r-^B.^i)-^nB} 



^2 



2c 



f-B]Bxr 



ihB] 



and, collecting terms, we get: 
dfl 

dt 4m 



dB 
'dt 



-^dB 
dt 



2m^ 



'x (S X n) - i^B 



(10) 



(11) 



The right hand side of Eq. ( pi] ) is not equal to ^B x fl = B x {r x IT) . The Bloch equations are replaced by 
operator equations whose consistent solution must treat the Heisenberg equation of motions for fl{t) as well as f(t) 
and n(t): 



dr ^ 

— — n/m 

dt ' 



~di 



an 

'di 



n,y(f) 



n,n /2m 



(12) 
(13) 



From these equations we see that energy can be absorbed by the system even under conditions of negligible or 
no coupling to the environment since dfl/ dt is no longer perpendicular to B. This power absorption in magnetic 
resonance systems containing electronic orbital angular momentum can also be interpreted as being due to the rate of 
change of the electric dipole moment in the presence of an electric field. This is evident from the relation: 




J\ ft _^ dB ^ d^E 
X r \ - U = -n ■ — = E ■ —— 
/ dt dt 



(14) 



Here the electron electric dipole moment is pE — qf — — er, its time derivative is dpE/dt = — eH/m, and the electric 
field is given by _E = —dA/dt. Clearly, additional power aborption may result due to coupling to a bath. 

The time derivative of the magnetic dipole moment in Eq. ( pT| ) can not be written simply in terms of the cross 
product of the external field and the dipole moment, hence a closed- form analytic solution of the operator equations 
in the Heisenberg representation is not available. Instead a perturbative solution can be used, as described below. 



IV. POWER ABSORPTION 

Substituting the magnetic field of Eq. (|^) into the Hamiltonian in Eq. (H), we find that the Hamiltonian can be 
separated into a sum of time-independent and time-dependent parts: 

H = Ho + Hi{t), (15) 

Ho = Vir) + 1^ - 5o^(ip. - VP.) + Bl^i.- + y^) , (16) 

2 

H,{t) = -Bi (e'"* + e-^'^*) h, ; h, = ^{yp, ~ zpy) + ^Bqxz , (17) 

2m 4m 
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where a term H2{t) of order B\ has been neglected. Realistic magnetic fields strengths impose a hierarchy between 
these terms because usually |i3o| ^ \Bi\- 

For comparison, two frequently used Hamiltonians which give the regular Bloch equation, Eq. are: — 
+ Hf{t), with: 



H'o = V{r) + BoAi 
2m 



(18) 

Hf{t)=H[^~Bi{e"^' + e-'^')h[; h[ ^ fi'^ , (19) 
obtained when all second order terms in the fields are neglected, as usual; and the Hamiltonian for a pure spin system: 

Hq = gsf^os- Bq, 



(20) 
(21) 



In both these cases for which the Bloch equation applies, the time- independent and time-dependent parts of the 
Hamiltonian commute: [H^,Hf] = 0, while [Ho,Hi] ^ 0. Either way, the time-independent Hamiltonian can be 
diagonalized in terms of a complete basis-set {|m)}: 



Ha\m) = E„,\m) 



(22) 



and the time-dependent electron wave function can be written (in the Schrodinger representation) as a linear combi- 
nation of these basis-set states: 



|*(t))=5]c„(t)^ 



-iEmt/h\ 



(23) 



with time-dependent coefficients, Cm{t), that satisfy a set of coupled equations, 

-"-"*(m|7?i|0Q , 



dcjYi i 



(24) 



where ujim = {Ei — Ejn)/fi- Power is absorbed by the system at the following rate: 

d 



V{u;) 



'dE' 




dt 


t 



dt 



EE 



E E c™^'^^ 



-iuJlmt 



{m\Hi\l) 



(25) 
(26) 

(27) 



Three different cases can be considered: 

• Pure spin system with no coupling to the environment. 

• Pure spin system with coupling to the environment. 

• Diamagnetic system (magnetic field induced magnetic moment due to orbital angular momentum) with no 
coupling to the environment. 

In the first case of a pure spin system with no coupling to the environment, as well as in the case of the approximate 
Hamiltonian leading to the regular Bloch equation; [H^,H^] — 0. As a result, the power absorption is identically 
zero. In this case it would be wrong to apply perturbation theory to the expansion of Eq. (|2^). This can be understood 
on very general terms: a time-dependent perturbation which commutes with the time-independent Hamiltonian can 
not induce transitions. 

Different magnetic systems (including in particular pure spin systems) with coupling to the environment were 
considered in the literature within two basically different approaches. In the first approach, Ti and T2 processes are 
included when considering the dynamics, and power absorption can thereby occur |l|J^. The Ti and T2 processes can 
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be described by relaxation terms added to the Bloch equation; the resulting dynamics can no longer be described by 
a Hamiltonian. 

In another approach introduced by Kubo and coworkers [Q, which we would like to compare to, it is assumed that 
a Hamiltonian completely defines the dynamics: 



H" = H'^ 



H[{t) 



7r-{yPz 
2m 



ZPy) 



(28) 
(29) 



Hq effectively incorporates the coupling to the environment; although the coupling is not explicitly specified, it is 
assumed that [Hq, H'l] ^ 0. The power absorption therefore does not vanish, and can be calculated using, for example, 
a linear- response theory Q| . 

We now show that the power absorption of a diamagnetic system with no coupling to the environment is very similar 
to power absorption of a pure spin system with coupling to the environment when the coupling to the environment 
is considered within linear- response theory as in Ref. Q . We show that in both cases the power absorption does not 
vanish and is given by a similar (though not identical) expression. The similarities and differences are derived and 
discussed below. The essential point is that as soon as [Hq,Hi] ^ 0, power absorption does not vanish, and, using 
the fact that for realistic magnetic fields i?i <C i?oj time-dependent perturbation theory can be applied. 

The first-order solution of Eq. (p3) with initial condition at t = 0, Cm(0) = 6™ is: 



Substituting this result into Eq. 
absorption is obtained: 



I 



1 



-i{u. 



1 



(30) 



^ — ^Im — UJrnl 

and applying the rotating wave approximation in the usual way, the power 



E„>E„ 



-|(n|/ii|m)|2(|6,„|2„|^^| 



sin (uj 



LO — LOn 



In the limit i — > oo, we find: 



(31) 



(32) 
(33) 



In the case of a pure spin system with coupling to the environment hi = jlx- Eq. ( p3| ) gives the static magnetic 
susceptibility and with Eq. ( p2| ) is identical to the results obtained by Kubo and Nagamiya Q in their Eqs. (18.4) and 
(19.17). We note that they have obtained this result within a different treatment, using a density matrix formalism, 
but under the same assumptions as those supporting the perturbative treatment we have used, namely, that the 
system is linear and in a steady-state. In the density matrix treatment, |6„ip and |6„p in Eq. (|3^) are replaced by 
the densities p{E,n) and p{En), respectively. 

We note here that this expression for the power absorption is correct only as long as a perturbative treatment is 
appropriate. As soon as a significant part of the population occupies the excited states a nonperturbative treatment 
would be necessary. The power absorption would necessarily cease as the system continues to absorb energy from the 
field and evolves to higher energy states if no coupling to the environment were present, but this is not contained in 
a perturbation theory or linear response approach. Here we study the power absorption only within the perturbative 
regime. We note that the time scale from turning on the fields till a significant change in the population could be 
noticed can in principle be quite large. 

In the case considered here, of no coupling to the environment but inclusion of orbital effects, hi — -^{yPz — zpy) + 

2 

-^Bqxz, while in the case of coupling to the environment but no orbital effects, hi — jl'^ = '^iVPz ~ zpy). Moreover, 
the spectrum and hence the resonance frequencies are different in these two cases. In either case, once the Bloch 
equation is modified, whether by relaxation terms or by terms of second order in the external fields, power absorption 
can ensue. Furthermore, while the numerical details of this power absorption may differ, the general behavior is the 
same in these two cases. The power absorption is given by Eq. ( ^2|) with a somewhat different quantum susceptibility in 
the two cases. In both cases, the quantum susceptibility is obtained by substituting the appropriate energy eigenvalues 
and coupling Hamiltonians into Eq. (^. 
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To summarize this section: focusing on the regime where a hnear response theory apphes, we have solved for the 
power absorption using time-dependent perturbation theory. We have derived Eqs. for the power absorption 

and for the susceptibihty. This result is similar to that obtained by Kubo and coworkers in Ref. Q but there are 
several differences: the coupling Hamiltonian hi is different, the quantum susceptibilities of the two systems are 
slightly different, and the spectral position is shifted. These are small corrections. Their order of magnitude is 
determined by the ratio between the magnetic energy and the total energy of the system. Hence, it seems that it 
would be difficult to experimentally observe a difference between the result for the susceptibility obtained using our 
approach and the standard result for cw applications where a linear response theory applies. The major difference 
is in justifying the linear response treatment. We have shown that it applies only if [Ho,Hi] ^ 0. This condition 
for a consistent application of perturbation theory holds in Kubo's derivation only if coupling to the environment is 
included, while for our case here no such coupling is required. 



V. POWER ABSORPTION OF AN ELECTRON BOUND IN A P-ATOMIC-ORBITAL 



As an example, we study power absorption of an electron bound by the Coulomb potential, V{r) = —Ze^ /{Attcqt), 
in a p-orbital {I = 1), and placed in the external time-dependent magnetic field of Eq. (||). Neglecting spin effects, 
the complete Hamiltonian can be written as above. As a basis set we choose the three wave functions: 



(f|V±> " ^ =^ fir) , (f|i/'o) = z f{r) 



where for principal quantum number rt = 2, 



z 



5/2 



— exp 



2ao 



(34) 



(35) 



and the Bohr radius is defined as ao = 47reo Tt? /{me^). These basis functions are eigenfunctions of the free Hamiltonian 
when no magnetic fields are present, as well as of -ffg, (but not of H or even Hq). The I = 1 manifold spanned by 
this basis set is not closed under the complete Hamiltonian. The orbital angular momentum / is not a good quantum 
number due to the terms in the Hamiltonian of second order in the external magnetic field, which do not commute 
with the angular momentum operator. The coupling to other manifolds is small, however, and the energy difference 
between the different manifolds is large. We therefore neglect the coupling to other manifolds, and consider the action 
of H only within the I = 1 subspace. A general solution of the time-dependent problem is then given by 

\^{t)) = co(t)e-^^°*/''|^o) +c+(t)e-^^+*/''|V'+) +c_(Oe-'^-*/''|V-) , 

and the time-dependent perturbative analysis of the coefficients of the previous section applies with m = 0, +, — . 
The diagonal matrix elements of H in this basis set are: 



12 



£ 



£ 15 

(^±|i7o|V±) = -g ±/ioSo + — (moSo)' 



where £ 



E± 

The level spacings are given by, 

hLU+ = E+ - Eq = fJ,oBo ^1 
hu!- = Eq — E^ = floBo ( 1 



>o , 

>0 , 



where hui- > assuming ShoBq < £. The oscillating field couples j-^o) and 1-0-1 



-Mo 



V2 



1 =F 12 



(36) 
(37) 

(38) 
(39) 

(40) 



The leading order contribution to the power absorption is: 



6 



X - \bo\') 6{c. - + (IfeoP - ^(c^ - . (41) 

The next order correction to this power absorption is: 

247r inoBo) , „ . . „ .2 
^i(^) = -^2 £ — (MoSo) (/io-Bi) 

X - |6oP) <5(^ - c._) + - l&oP) ^(^ - ^+)} . (42) 

Higher order terms are of third order in Bi. At temperature T — 0, bo = = and 6_ = 1, and we obtain: 



7^ (mo^o) ifJ'oBi f 
n 





w - 


Mo-Bo A 








a V 


f )\ 



(43) 



Note that for reahstic experimental conditions, (/io-Bi) is small enough to justify the perturbative treatment and as 
a result V{uj) is very small. 

Realizable magnetic fields are, for example, i?o = 5 T and Bi = 5 x 10"'* T. The relative strength of the correction 
term Vi as compared to the leading order term Vq, is determined by the dimensionless parameter 24/ioi?o/£, which 
for Bq = 5 T and Z — 6 (carbon) is equal to 1 x 10~''^. Vi as well as the higher order terms are thus but small 
corrections to VoicLi). 



VI. CONCLUSIONS AND SUMMARY 



In this paper we have studied the dynamics of the gauge-invariant operator for the magnetic dipole moment /i of 
an electronic system in an external magnetic field B. In addition to the intrinsic magnetic moment of the electron it 
includes a field induced magnetic moment arising from the electronic orbital angular momentum, i.e., a diamagnetic 
contribution. We have derived the Heisenberg equation of motion for this magnetic moment in a general magnetic field 
and shown that this modified Bloch equation is coupled to two other Heisenberg equations for the position and velocity 
operators. While we could not solve this set of equations exactly and in closed form, we have solved it perturbatively 
using time-dependent perturbation theory and calculated the power absorption of the magnetic moment in a time- 
dependent magnetic field under magnetic-resonance conditions. Under the assumption of a steady-state, we have 
found an expression for the quantum susceptibility in terms of matrix elements of a coupling Hamiltonian hi. Finally, 
we have explicitly calculated the integrated power absorption for the example of a p-orbital electron and found that 
it is given by Eq. (^) . 

Several generalizations can be easily considered. Upon generalizing to include electron spin, the Hamiltonian and 
the magnetic moment operator become, 

H = V{r) + ^-ll,-B , (44) 



/7 = — [r x'R + gsTis) = — r x H - .gsMos , 
2m 2m 



(45) 



respectively, where fig — —gsi-iQS, and gs is the electron spin gyromagnetic ratio (g^ ~ 2.002). The last term on the 
right hand side of Eq. (BS) is the intrinsic spin magnetic moment operator of the electron, fl^. Incorporation of this 



term into the modified Bloch equation yields the following term on the right hand side of Eq. ( pT| ): 
i.e., the generalized Bloch equation becomes. 



dfl 



e 

4m 



dB 
'dt 



,dB_ 
' dt 



2m? 



'x (B X H) - ifiB 



h 



■B X s 



-B X 



(46) 



Inclusion of spin-orbit coupling of the form a{r)s-l into the Hamiltonian is also straightforward. The modified Bloch 
equation then has additional terms on the right hand side proportional to a{r)sx I and a{r){s- r){B ■ r). The latter 
term is due to the magnetically induced magnetic moment. Furthermore, one can also readily include ligand field 
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effects, exchange coupling effects, etc., within the formulation considered here. Ligand field terms in the Hamiltonian 
commute with the magnetically induced magnetic moment, as do spin exchange terms, but orbital angular momentum 
exchange coupling terms do not, and therefore contribute additional terms to the extended Bloch equation due of the 
magnetically induced magnetic moment. 

How do our results differ from previous treatments of magnetic resonance phenomena? It is useful to divide such 
a comparison into two subjects: the equations of motion and the mechanism of power absorption. 

We have shown that the Bloch equation, Eq. (^, is not the correct Heisenberg equation for a magnetic moment /2 
in the magnetic field B when orbital angular momentum is present. A modified Bloch equation, Eq. (|l^), is obtained. 
Upon inclusion of spin, the factor of gs in the term —ii^B x fl^, which is equivalent to the last term on the right 
hand side of Eq. (^) , further precludes writing a Bloch equation of the form of Eq. (^) with jl replaced with jli + jl^. 
The modified equation can not be written in terms of fl and B alone. It contains the position and kinetic velocity 
operators and therefore must be solved together with dynamical equations for these operators. In general, dfl/dt is 
not perpendicular to B, and fl oscillates out of phase with B. This results in a non-vanishing power absorption. 

In the usual picture, correct for pure spin systems, power is absorbed only because of the coupling of the spin 
system to the environment. Without this coupling the magnetic moment fl would oscillate in phase with B, dfl/dt 
would be perpendicular to B at all times, and no power would be absorbed. Coupling to the environment dephases 
the oscillations of fl and B and induces a nonvanishing component of dfl/dt along B. In the case considered here, of 
a magnetic moment that explicitly contains a diamagnetic contribution, the power absorption does not vanish even 
when coupling to the environment is not included because the magnetic moment has a component which oscillates out 
of phase with the external field as a result of non-linear terms in the interaction Hamiltonian arising from a magnetic 
field induced magnetic moment. The phase difference between oscillations of the magnetic field and oscillations of 
the magnetic dipole moment of a bound electron in this field, is an intrinsic feature of the system and need not be 
induced by additional couplings to a bath that gives rise to dissipation. This power absorption can be interpreted, as 
evident from Eq. ([ij) , as being due to the rate of change of the electric dipole moment in the presence of an electric 
field. Additional power aborption may result due to coupling to a bath. 

For cw power absorption, dissipation of energy from the atomic system to the environment is required in order to 
ensure that the system stays in thermal equilibrium, so that power absorption can be maintained under steady-state 
conditions. Hence, for truly cw power absorption, as opposed to pulsed magnetic resonance experiments, the two 
systems depend on coupling to the environment to maintain a steady state and eventually will behave similarly. The 
results presented here indicate that diamagnetic systems would behave differently from pure-spin systems during the 
initial time when the turning on of the fields occur. For pure spin systems the time scale for initiation of power 
absorption would be set by the coupling to the environment while for diamagnetic systems it would be an intrinsic 
time scale indifferent to the environment. Without a relaxation mechanism, no power would be absorbed by a pure 
spin system, while the diamagnetic system will absorb power upon the turning on of the fields only until the population 
of levels reach a new equilibrium. We note that as soon as the population is changed significantly, second and higher 
order effects become important and a first-order time-dependent perturbation theory would not suffice. 
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